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Beschreibung 

Anwendungsgebiet 

Die mechanischen Eigenschaften von biologischem Ge- 
webe (z. B. Elastizitatsparameter) sind fur die Beurteilung 
des Zustandes des Gewebes von groBem Interesse. In der 
medizinischen Diagnostik deuten Veranderungen der Elasti- 
zitatseigenschaften auf histoLogische und u. U. pathologi- 
sche Veranderungen hin. Allgemein bekannt sind Prozesse 
wie die Bildung von Geschwulsten und Verhartungen 
("Knoten"), die z. B. manuell tastbar sind ("Palpation"). In 
der Landwirtschaft ist. fur die Beurteilung der Qualitat des 
Fleisches von Schlachtvieh die Kenntnis mechanischer Ge- 
webeeigenschaften ebenfalls von Interesse. 

Stand der Technik 

Die sogenannte Tastbefundung ist ungenau und unemp- 
findlich. Wesentlich besser ist in dieser Beziehung die soge- 
nannte Elastographie, bei der elastische Gewebeeigenschaf- 
ten technisch erfaBt und z. B. in Form von Schnittbildem 
qualitativ oder quantitativ visualisiert werden. Dabei be- 
dient man sich hauptsachlich des Ultraschalls, wie er als 
bildgebendes Verfahren in der medizinischen Diagnostik an- 
gewandt wird, allerdings in einer modifizierten Betriebsart. 
In zeitlich nacheinander aufgenommenen Ultraschallbildern 
konnen geringste Verschiebungen oder Verfcrmungen in- 
nerhalb der dargestellten Gewebes truktur durch Auswertung 
der Bildsequenzen erfaBt und ausgewertet werden. Wird auf 
einen Gewebebereich ein mechanischer Druck ausgeubt, der 
eine Verformung des Gewebes zur Folge hat, so verformen 
sich Bereiche mit unterschiedlichen Elastizitatseigenschaf- 
ten verschiedenartig. Das Elastographiesystem wertet diese 
Verformungen durch den numerischen Vergleich der Einzel- 
bilder aus und stellt die unterschiedliche Dehnungen im Bild 
dar. Die notwendige Kompression des Gewebes wird extern 
durch den Ultraschallwandler provoziert oder intern durch 
Atmung oder Puis. Die Kompression ist nur gering und be- 
tragt bei Anwendung des ublichen diagnostischen Ultra- 
schalls nur Bruchteile von Millimetern, Wichtig dabei ist 
eine quantitative Kontrolle des Kompressionsdruckes. Der 
Druck wird in Richtung der Schallausbreitung provoziert. 

Ein Verfahren der Ultraschall-Elastographie von Korper- 
gewebe ist. erstnialig in einem Aufsatz von J. Ophir et al iin 
Jahrc 1991 [1], [2] bcschricbcn wordcn. Dabei werden Ul- 
traschalibilder bzw. die korrespondierenden hochfrequenten 
Ultraschallechosignale, aus denen die ublichen Ultraschall- 
biider in dern Gerat entstehen, so ausgewertet, daB Verschie- 
bungen des Korpergewebes zwischen zwei, mit verschiede- 
ner Kompression aufgenommennen Gewebebildern berech- 
net werden. Auf diese Weise lassen sich die bereits erlauter- 
ten Ruckschlusse auf die Elastizitat des Organs bis hin zu ei- 
ner quantirativen Abbildung des Elastizitatsmoduls erzielen. 

Die hochfrequenten Echosignale aus einem komprimier- 
ten Gewebeabschnitt erreichen den Ultraschallwandler zu 
einem Zeitpunkt der von dem MaB der Kompression des be- 
treffenden komprimierten Gewebeabschnitts abhangt. Die 
Hauptaufgabe in der Auswertung der hochfrequenten Ultra- 
schallechosignale zur Berechnung der Dehnung besteht da- 
her darin, Zeitverschiebungen in kurzen zeitlichen Interval- 
len der hochfrequenten Ultraschalldaten zu berechnen. 
Diese Zeitverschiebungen zwischen korrespondierenden 
Echosignalen werden in alien durch das Ultraschallsystem 
abgebildeten Gewebebereichen in Form eines zweidimen- 
sionalen Bildes berechnet. Fur die Berechnung einer Zeit- 
verschiebung wird die Kreuzkorrelationsfunktion zwischen 
den hochfrequenten Echodaten verwendet. Einerseits kann 



die Zeitverschiebung durch Sue hen des Maximums der 
Kreuzkorrelationsfunktion ermittelt werden [1]. Anderer- 
seits kann die Zeitverschiebung aus der Phase der Kreuzkor- 
relationsfunktion an irgendeiner S telle, i. d. Regel an der 
5 Steile Null, geschatzt werden [3], [4]. 

Nach der Berechnung der Zeitverschiebung kann die lo- 
kale Dehnung durch Gradientenbildung oder mit. Hilfe ein- 
facher linearer Filter [6] berechnet und abgebildet werden. 
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Nachtelle des Standes der Technik 



Fur alle Verfahren zur Berechnung der Zeitverschiebun- 
gen wird die Kreuzkorrelationsfunktion korrespondierender 
Intervalle der Echosignale aus dem gleichen Gewebeab- 
15 schnitt unter verschiedener Kompression verwendet. Die 
Auswertung der Kreuzkorrelationsfunktion ist aufgrund der 
notwendigen Integration Uber das ganze Intervall zeitauf- 
wendig. Die bisher vorgeschlagenen Verfahren lassen sich 
in zwei Gruppen unterteiien: 

20 

1) Kreuzkorrelationsmethoden 

Verfahren, die das Maximum der Kreuzkorrelationsfunk- 
tion durch eine vollstandige Suche oder Einschachtlung er- 

25 mitteln. Dazu miissen die Echosignale in der Regel mit einer 
sehr hohen Abtastrate abgetastet. werden, um auch kleine 
Zeitverschiebungen genau messen zu konnen. Diese Verfah- 
ren sind auBerst zeitaufwendig, so daB sie nicht in ein Echt- 
zeitsystem oder Online-System implementiert werden kon- 

30 nen. 

2) Phasenbasierte Verfahren 

Verfaliren, die Zeitverschiebungen aus der Phase der 
35 Kreuzkorrelationsfunktion der komplexwertigen Ultra- 
schallechosignale schatzen. Nachteile dieser Verfaliren sind 
mogliche Ungenauigkeiten und das Auftreten von Mehrdeu- 
tigkeiten bei groBen Zeitverschiebungen. Bisher konnten 
diese Mehrdeudgkeiten nur durch zusatzliche zeitin tensive 
40 zweidimensionale V^r- oder Nachverarbeitungsschritte ge- 
wahrleistet werden. 

Aufgrund der Nachteile der bisherigen Verfahren ist es 
nicht moglich, Ultraschallechodaten hinreichend genau in 
Echzeit zu verarbeiten und somit Dehnungsbilder in Echzeit 
45 darzustellen. In der Vergangenheit hat sich aber gezeigt, daB 
bei Ultraschallabbildungsvcrfahrcn wie dem klassischcn B- 
Bild oder der Dopplersonographie die Echzeitfahigkeit der 
Systeme besonders wichtig fur eine breite Nutzung dieser 
Systeme ist. 

50 Ein weiteres Problem ist, daB die Veranderungen der Si- 
gnallaufzeiten (Puls-Echo-Sequenz) die sich durch die 
Kompression des Gewebes ergeben, von der Entfernung 
zwischen dem Gewebeareal, welche die Echos hervorrufen, 
und dem Ultraschallwandler abhangen, d. h. nicht konstant 

55 sind und mit zunehmender Entfernung monoton ansteigen. 
Durch Schwankungen der Einhullenden der Echosignale 
kommt es dalier zu Ungenauigkeiten in der Verschiebungs- 
berechnung, da diese ein Intervalls bestimmter Lange be- 
nutzt. Diese wurden in der Vergangenheit u. a. durch Log- 

60 arithmieren der hochfrequenten reellen Echosignale [7] ver- 
mindert. Nachteilig wirkt sich hierbei aber die Veranderung 
der Phasen der Signale aus. 



65 



Aufgabe der Erfindung 

Aufgabe der Erfindung ist es daher, eine Anweisung zum 
Handeln zu entwickeln, mit dem Ultraschallechodaten hin- 
reichend genau moglichst schnell (z. B. in Echtzeit) ausge- 
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wertet werden konnen. Das Verfahren sollte die gleiehe Ge- 
nauigkeit wie die Kreuzkorrelationsverfaliren durch eine 
moglichst recheneffiziente Signalverarbeitung erreichen, 

Losung der Aufgabe 

In phasenbasierten Verfahren wird die Phasenverschie- 
bung zwischen zwei Echosignalen ais MaB far die Zeitver- 
schiebung benutzt. Fur monofrequente zeitverschobene Si- 
gnale ist die Phase der Kreuzkorrelationsfunktion eine li- 
neare Funktion der zeitlichen Echoverschiebung. Bei der 
tatsachlichen, d. h. gesuchten Zeitverschiebung hat die 
Phase der Kreuzkorrelationsfunktion eine Nullstelle. Die 
Sleigung entspricht der Kreisfrequenz der monofrequenten 
Signale. Dieser Zusammenhang kann fur die Berechnung 
der Zeitverschiebung aus der Phase der komplexen Kreuz- 
korrelationsfunktion an einer beliebigen Stelle genutzt wer- 
den. Die Phase der Kreuzkorrelationsfunktion kann aller- 
dings aufgrund der Mehrdeutigkeit der Arcustangensfunk- 
tion numerisch nicht eindeutig ermittelt werden. AuBerdem 
werden in Ultraschall-Puls-Echosystemen moglichst breit- 
bandige Signale verwendet. Fur bandbegrenzte Signale gilt 
obiger Zusammenhang nur approximative Diese Tatsachen 
fuhren zu den genannten Problemen der Ungenauigkeit und 
der Mehrdeutigkeit phasenbasierter Verfahren. 

Beide Probleme (Mehrdeutigkeit, Ungenauigkeit bei 
Breitbandsignalen) konnen gelost werden, indem iterativ die 
Phase der Kreuzkorrelationsfunktion fur die Suche des Ma- 
ximums der Kreuzkorrelationsfunktion verwendet wird. Die 
vom Ultraschallwandler aufgenommenen Echodaten sind 
gemafi Abb* 1 in A-Linien organisiert. Die Berechnung von 
Zeitverschiebungen erfolgt an verschiedenen diskreten 
Punkten im Gewebe. In der Elastographie ist die Zeitver- 
schiebung, die an diesen diskreten Punkten berechnet wer- 
den konnen ein stetige, monotone Funktion uber der Hefe. 

Die Zeitverschiebungen werden sukzessiv fur jede A-Li- 
nie beginnend vom dem Ultraschallwandler nachstgelegen- 
den Punkt bis zum tiefsten Punkt berechnet. Fur die Berech- 
nung der Zeitverschiebungen der Echos eines Punktes wird 
die folgende Iteration durchlaufen: 

1. ) Die Iteration wird begonnen bei der geschatzten 
Zeitverschiebung %, die fur den in der gleichen A-Linie 
dem Ultraschallwandler naherliegenden Punkt ermit- 
telt wurde. Fur die Berechnung der dem Ultraschall- 
wandler nachstgclcgcncn Punkte wird die Iteration bei 
Null begonnen. 

2. ) Die Phase der Kreuzkorrelationsfunktion der kor- 
respondierenden Ultraschallechodaten an der geschatz- 
ten Zeitverschiebung % wird ermittelt. Die ermitteite 
Phase ist aufgrund obiger Eriauterungen ungefahr pro- 
portional zum Abstand der geschatzten Zeitverschie- 
bung zur wirklichen Zeitverschiebung. Die Proportio- 
nalitatskonstante entspricht in etwa der Schwerpunkts- 
frequenz des verwendeten Ultraschallwandlers (Kreis- 
frequenz). Die Schatzung der Zeitverschiebung X kann 
daher verbessert. werden, indem der verbleibende Pha- 
senunterschied durch Teilen durch die Schwerpunkts- 
frequenz des Ultraschallwandlers in eine Zeitverschie- 
bung umgerechnet und die bisherige Schatzung der 
Zeitverschiebung entsprechend verandert wird. 

Durch Punkt 1.) wird gewahrleistet, daB die Kreuzkorre- 
lationsfunktion stets in der Nahe der wirklichen Zeitver- 
schiebung ausge wertet wird. So kann gewahrleistet werden, 
daB verbleibende Phasenunterschiede eindeutig in Zeitver- 
schiebungen umgerechnet werden konnen. 

Durch eine mehrmaliges iteratives Wiederholen von 2.) 



wird die Ungenauigkeit bei der Umrechnung von Phasenun- 
terschieden in Zeitverschiebungen bei der Verwendung 
breitbandiger Signale iterativ minimiert. Insgesamt muB mit 
diesem Verfahren die Zeitverschiebung die Kreuzkorrelati- 

5 onsfunktion lediglich wenige Male (weniger als 10 Mai) fur 
die hinreichend genaue Schatzung einer Zeitverschiebung 
ausgewertet werden. Das vorgestellte Verfahren ist somit in 
zwei Beziehungen iterativ: Zeitverschiebungen werden 
durch wiederholtes Verwenden von berechneten Zeitver- 

10 schiebungen gleicher oder vorangehender Bildpunkte be- 
rechnet. Da in der Elastographie auch sehr kleine Zeitver- 
schiebungen geschatzt werden miissen, die in der Regeln 
unterhalb der Abtastzeit liegen, miissen die Echosignale ge- 
eignet interpoliert werden. Bei dem vorliegenden Verfahren 

15 werden wie in [3] und [4] fur die Berechnung der Kreuzkor- 
relationsfunktion daher die komplexen, heruntergernischten 
Basisbandsignale verwendet. Da diese aufgrund der Band- 
paBbegrenzung von Ultraschall-Echosignalen eine deutlich 
geringere obere Grenzfrequenz aufweisen, reicht hier eine 

20 recheneffiziente lineare Interpolation aus [5]. 

Ungenauigkeiten der geschatzten Phasenverschiebungen 
werden in diesem Verfahren dadurch vermindert, daB die 
Amplituden der komplexen Basisbandsignale logarithmiert 
werden, nachdem sie zuvor um eins erhoht wurden. Die 

25 Phasen werden nicht verandert. Dadurch kann die logarith- 
mische Kompression auch fur phasenbasierte Verfahren be- 
nutzt werden. 

Vorteite der Erfindung 

30 

Das Verf aliren erreicht im Gegensatz zu bisherigen pha- 
senbasierten Verfahren die gleiehe Genauigkeit wie Kreuz- 
korrelationsmethoden. Dazu muB die Kreuzkorrelations- 
funktion nur einige Male (weniger als 10 Mai) ausgewertet 

35 werden. Auf diese Weise ist eine Echtzeitauswertung der 
Echosignale zur Berechnung der Gewebedehnung bereits 
jetzt mit geeigneten digitalen Signalverarbeitungsanlagen, 
wie sie in modernen digitalen Ultraschallgeraten zur Verfu- 
gung stehen oder separat als kommerziell erhaltiich sind, 

40 moglich. 
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Zeichnungen dargestellt und werden im folgenden naher be- 
schreiben: 

Zu Abb, 2) 

5 

Schematische Darstellung eines System zur Aufnahme, 
schnellen Berechnung und von Darstellung von Dehnungs- 
bildern 

Das System zur Aufnahme, schnellen Berechnung und 10 
Darstellung von Dehnungsbildern besteht aus drei Funkti- 
onseinheiten: 

a) ein handelsubliches Ultraschallsystem mit verschie- 
denen Sonden. Diese konnen Ultraschallechodaten mit 15 
verschiedenen Scanprinzipien aufnehmen (Sektor- 
Scan, Parallel-Scan). Ebenso spielt die Art der Fokus- 
sierung (elektronisch oder durch einen Einzelwandler 
mit konstanteni Fokus) keine Rolle. Die hochfrequen- 
ten Echodaten miissen in digitalisierter Form vorlie- 20 
gen. Bei Systemen mit analogem Beamformer, oder bei 
Systemen mit Einzelwandler miissen die analogen 
hochfrequenten Echodaten direkt nach dem Beamfor- 
mer, bzw. nach dem Einzel wander abgegriffen und di- 
gitalisiert werden. Bei Systemen mit digitalem Beam- 25 
former konnen die digitaien Echodaten direkt abgegrif- 
fen werden. 

b) eine Kecheneinheit zur schnellen Berechnung der 
Dehnung. Dies ist das Kernstuck der Erfindung. Diese 
Recheneinheit kann sich innerhalb oder auBerhalb des 30 
Ultraschailgerates befinden. Die Recheneinheit kann 
durch eine digitale Schaltung oder aber durch eine 
DSP-Karte, einen anderen Rechner, eine Recheneinheit. 
des Ultraschailgerates realisiert werden (siehe Abb. 3) 

c) eine Darstellungseinheit Diese Darsteliungseinheit 35 
kann sich ebenfalls innerhalb oder auBerhalb des Ultra- 
schailgerates befinden. Es konnen die Darsteliungs- 
moglichkeiten des Ultraschailgerates verwendet wer- 
den. Bei AnschluB eines externen Rechner als Rechen- 
einheit konnen grafische Darstellungsmoglichkeiten 40 
dieses Rechners verwendet werden. 

Die Berechnung der Dehnungsbilder kann in Echtzeit, 
online oder offline nach der Aufnahme der Ultraschallcho- 
daten erfolgen. Fur die Darstellung der Dehnungsbilder 45 
kann das Dchnungsbild komplctt dargestellt werden, ein 
Teil des B-Bildes in einem bestimmten Bereich (Region of 
Interest) durch das Dehnungsbilder iiberblendet oder einge- 
farbt werden. 

50 

zu Abb. 3 

Realisierung der Recheneinheit 

Die Recheneinheit kann in folgender Weise z. B. als Digi- 55 
talschaltung realisiert werden: zunachst miissen die hochfre- 
quenten Echodaten quadraturdemoduliert werden (1), Die 
resultierenden Daten sind abwartsgemischt. Sie werden Ba- 
sisbanddaten oder IQ-Daten genannt. Die Amplituden der 
IQ-Daten werden nun wahlweise logarithimsch kompri- 60 
miert (2). 

Die entsprechend aufbereiteten komplexen Echosignale 
(IQ-Daten) werden in einem Zwischenspeicher (3) abge- 
speichert. Dieser Zwischenspeicher muB die IQ-Daten je- 
weils zweier zeitlich aufeinanderfolgender B-Bilder spei- 65 
chem konnen, Er kann daher nach dem FTFO-Prinzip im- 
mer, wenn die Echodaten eines neuen Bildes am Eingang 
anliegen, die Echodaten des vorletzten Bildes loschen. 



Die Berechnung der Zeitverschiebung erfolgt in einer ge- 
sonderten Einheit (11). Diese kann in einer digitaien Schal- 
tung mehrfach vorhanden sein. Dadurch konnen die Zeitver- 
schiebungen mehrerer A-Linien gleichzeitig berechnet wer- 
den. Der Ablauf wird von einer globalen Ablaufsteuerung 
(12) gesteuert. Diese Ablaufsteuerung steuert die Berech- 
nung der Zeitverschiebung an verschiedenen diskreten Posi- 
tioned Die Zeitverschiebungen werden gemaB Abb. 1 an 
verschiedenen diskreten Bildpunkten berechnet. Verschie- 
dene A-Linien konnen parallel oder nacheinander abgear- 
beitet werden. Bei der Berechnung der Zeitverschiebungen 
einer A-Linie wird mit dem den Wandler am nachsten lie- 
genden Punkt begonnen und danach sukzessive Zeitver- 
schiebungen an tiefer liegenden Punkten berechnet. Dieser 
Ablauf muB von der globalen Ablaufsteuerung (12) kontrol- 
liert werden. 

Die Einheit zur Berechnung einer Zeitverschiebung (11) 
enthalt folgende Elemente: 

Die komplexe Kreuzkorrelationsfunktion wird an der durch 
den Integrator (9) vorgegebenen Stelle ausgewertet (5). Der 
Integrator (9) dient als Speicher fur die geschatzte Zeitver- 
schiebung wahrend der Iteration, Sein Startwert muB von 
der globalen Ablaufsteuerung gesetzt werden. Der Startwert 
fur die dem Ultraschallwandler nachstgelegenen Punkte ist 
Null. Der Startwert fur alle anderen Punkte ist die Zeitver- 
schiebung, die fur den dem Ultraschallwandler nahergelege- 
nen Punkt berechnet wurde. Die Lage des Intervalls der Da- 
ten, die fur die Berechnung der Kreuzkorrelationsfunktion 
verwendet wird, wird von der globalen Ablaufsteuerung 
(12) bestimmt. Dieses ist ein Intervall konstanter Lange um 
den Punkt an dem die Zeitverschiebung berechnet werden 
soil. 

Damit die Kreuzkorrelationsfunktion auch an Stellen be- 
rechnet werden kann, die nicht ganzzahlige Vielfache der 
Abt astzeit sind, miissen die Daten linear von einem Interpo- 
lator (5) interpoliert werden. Dieser greift auf die Basis- 
band-Echodaten der gleichen A-Linie der beiden B-Bilder 
im Bildspeicher unter unterschiediicher Kompression zu. In 
(5) wird somit die Kreuzkorrelationsfunktion der Basis- 
banddaten berechnet. Um die Kreuzkorrelationsfunktion der 
komplexen hochfrequente Daten zu erhalten, muB die Phase 
des Ergebnisses um das Produkt der Kreisfrequenz des Mi- 
schers (1) mit dem Zeitpunkt an dem die Kreuzkorrelations- 
funktion ausgewertet wird (9) in negative Richtung gedreht 
werden (6). Die Phase dieser komplexen Zahl wird nun be- 
rechnet (7) und mit dem rcziprokcn Wert der Schwcrpunkts- 
frequenz des Wandlers (Kreisfrequenz) multipliziert (8). 
Das Ergebnis wird zum Wert des Interpolators (9) hinzuad- 
diert Der Ablauf zwischen (5) und (9) muB einige Male 
(1-10) wiederholt werden. Dies wird von einer Ablaufsteue- 
rung (10) kontrolliert Das Ergebnis der Iteration wird an- 
schlieBend in einem Ergebnisbildspeicher (13) abgelegt. Die 
Stelle, an der das Ergebnis dort abgelegt wird, wird von der 
globalen Ablaufsteuerung kontrolliert.. Der Inhalt des Bild- 
speichers kann anschlieBend zur Anzeige gebracht werden. 

Samtliche Elemente dieser Recheneinheit konnen auf ver- 
schiedene Weise technisch realisiert werden: 

1 . Jedes oder einige Element kann durch eine digitale 
Logikschaltung realisiert werden 

2. Jedes oder einige Elemente konnen von einer oder 
mehreren DSP-Schaltungen realisiert werden. 

3. Jedes oder einige Elemente konnen von einem 
Rechner innerhalb oder auBerhalb des Ultraschailgera- 
tes realisiert werden. 
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Patentanspruche 

1 . Ein Verfahren und ein System zur Berechnung der 
Dehnung elastischen Gewebes mit diagnostischem Ul- 
traschall, dadurch gekennzeich net daB die Zeitver- 5 
schiebung durch iterative Auswertung der Phasen der 
Kreuzkorrelationsfunktion zugehoriger EchosignaLe 
geschatzt werden. 

2. Ein System naeh Ansprucb 1, dadurch gekennzeich- 
net, daB die iterative Auswertung der Phasen der 10 
Kreuzkorrelationsfunktion an transducernahen Gewe- 
bebereichen beginnt, bei denen die Nullstellensuche 
eindeutig ist und, daB hier gewonnene Ergebnisse bei 
der Berechnung von Verschiebungen bzw. Phasennull- 
durchgange in weiter entfernten Gewebebereichen in 15 
der Weise vorteilhaft genutzt werden konnen, daB die 
Kenntnis der vorher gewonnenen Phaseninforniationen 
bei entfernten Bereichen vorhandene Uneindeutigkei- 
ten beseitigt. 

3. Ein System nach An spruch 1, dadurch gekennzeich- 20 
net, daB die notwendigen Interpolalionen hochfrequen- 
ter Echosignale mit Hilfe einer linearen Interpolation 
von abwartsgemischten Echosignalen (Basisbandsi- 
gnalen) erfolgL 

4. Ein System nach Anspruch 1, dadurch gekennzeich- 25 
net, daB vor der die Genauigkeit der Berechnung von 
Zeitverschiebungen dadurch verbessert wird, daB die 
Amplitude der komplexen Echosignale (unter Beibe- 
baltung der Phase) logarithmisch komprimiert werden. 

30 

Hierzu 3 Seite(n) Zeichnungen 
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